The Cole-Cole Law for Critical Dynamics in Glass-Forming Liquids 
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Within the mode-coupling theory (MCT) for glassy dynamics, the asymptotic low-frequency expansions for 
the dynamical susceptibilities at critical points are compared to the expansions for the dynamic moduli; this 
shows that the convergence properties of the two expansions can be quite different. In some parameter regions, 
the leading-order expansion formula for the modulus describes the solutions of the MCT equations of motion 
outside the transient regime successfully; at the same time, the leading- and next-to-leading order expansion 
formulas for the susceptibility fail. In these cases, one can derive a Cole-Cole law for the susceptibilities; and 
this law accounts for the dynamics for frequencies below the band of microscopic excitations and above the 
high-frequency part of the a-peak. It is shown that this scenario explains the optical-Kerr-effect data measured 
for salol and benzophenone (BZP). For BZP it is inferred that the depolarized light-scattering spectra exhibit a 
wing for the a-peak within the Gigahertz band. This wing results from the crossover of the von Schweidler- 
law part of the a-peak to the high-frequency part of the Cole-Cole peak; and this crossover can be described 
quantitatively by the leading-order formulas of MCT for the modulus. 



I. INTRODUCTION 

During the past 15 years, several new spectrometers have 
been introduced for the study of the glassy dynamics of liq- 
uids. The evolution of this complex slow dynamics upon de- 
creasing the temperature T or increasing the density p has 
been documented for many systems for times t, which ex- 
ceed the natural time scale f m ; c for condensed-matter motions 
by three or more orders of magnitude. Similar progress has 
been made for molecular-dynamics simulations of liquid mod- 
els. In parallel to these experimental activities, a theory for 
the evolution of glassy dynamics has been developed which 
is referred to as mode-coupling theory (MCT). This theory 
is based on regular equations of motion for a set of auto- 
correlation functions. The MCT equations lead to fold bifur- 
cations for the correlators' long-time limits if some control 
parameter like T reaches a critical value T c ; this bifurcation 
describes a transition from a liquid to an amorphous solid. 
The distance of the control parameter from the critical value, 
say e = (T c — T)/T c , can be used for the discussion of the 
bifurcation dynamics as a small parameter. For £ tending to 
zero and times increasing to infinity, it is possible to calcu- 
late asymptotic solutions of the MCT equations. The leading- 
order results provide a set of general formulas, which explain 
the qualitative features of the bifurcation scenario. Many fits 
of data with these general formulas have been studied in or- 
der to test the relevance of the MCT for the explanation of the 
experimental facts Q|. 

A set of general MCT results, which is of main interest in 
this paper, concerns the critical dynamics. This dynamics is 
described by auto-correlation functions (j)(f) for control pa- 
rameters at the bifurcation point, say T = T c . Equivalently, 
one can consider the corresponding loss spectra %"(co). These 
are the products of the frequency co and the Fourier-cosine 
transform of §(t). The central asymptotic formula is speci- 
fied by (1) a positive number f c , which is called the plateau, 
(2) a positive amplitude, say A, and (3) the critical exponent a, 
obeying < a < 0.396, lim,_ >00 f fl [(|)(f ) - f] = A. This formula 
is equivalent to lim ffl ^o X'^ 10 )/ 10 " = sin(7ta/2)r(l — a)A, with 
r denoting the gamma function. For a given transition point, 



all correlators are specified by the same exponent a, but differ- 
ent critical points can differ in their value for a. The leading- 
order long-time result for the correlator describes a power- 
law decay: (j)(f) — f c °< l/t". Equivalently, the leading-order 
result for the low-frequency critical loss spectrum is given 
by a power law variation %"(co) « co". For states near the 
transition point, (j)(f) and %"(&) can be replaced by their re- 
spective critical functions for short times, t -C f a , or large 
frequencies, C0f o 3> 1. The time scale f a is the same for all 
correlators and diverges for states approaching the transition 
point. For t f a and cof a ^ 1, the correlators and spectra de- 
pend sensitively on e; for shorter times, t <C t a and C0f o ^> 1, 
the correlators and spectra depend on e smoothly. An co" 
spectrum was identified first for the glass-forming molten salt 
0ACa(NO3) 2 0.6K(NO3) (CKN) in data obtained by neutron- 
scattering spectroscopy \2]. This system was also used to 
document for the first time the evolution of glassy dynam- 
ics within the full Giga-Hertz band |3]: Using depolarized 
light-scattering spectroscopy, a spectrum compatible with the 
co" law was found extending from 1GHz to 400GHz. The 
t~ a decay in the time domain was measured first for density 
correlators by photon-correlation spectroscopy for a colloidal 
suspension of hard spheres |4] with the density as a control 
parameter. 

Other glass-forming systems that can be studied experi- 
mentally are found in many van-der Waals liquids. The nat- 
ural time scale for inter-molecular vibrations is a picosec- 
ond. For normal liquid behavior, one expects correlations 
to decay to zero for times around some picoseconds. The 
normal-liquid excitation spectra extend from, say, 0.5THz to, 
say, 5THz. Glassy-dynamics spectra for several systems have 
been measured by depolarized light-scattering spectroscopy, 
and the data were shown to be consistent with the MCT bifur- 
cation scenario. For example, spectra for toluene have been 
fitted successfully with leading-order asymptotic results for 
frequencies between 0.5GHz and ITHz and for temperatures 
decreasing from T = T c + 140K to T c + 10K; T c w 150K H. 
For all these systems, the weight of the glassy-dynamics part 
of the spectra is large compared to the part of normal-liquid 
dynamics, i.e., the so-called a-peak is large. In agreement 
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with the MCT prediction for such situations, the amplitude A 
for the co a spectrum is small, and the low-frequency contribu- 
tions of the normal-liquid dynamics affect the region of the 
expected co" behavior of the spectra. As a result, no frequency 
interval has ever been identified for a van-der-Waals liquid, 
where an of spectrum can be identified explicitly. 

Torre et. al |6] have introduced optical-Kerr-effect (OKE) 
spectroscopy as a technique for the study of glassy dynam- 
ics. The measurement provides the response function %(t) oc 
— 3f0(f) for the same probing variable, which is studied in 
depolarized light-scattering experiments. The Fourier-sine 
transforms of %(t) are proportional to the loss spectra %"(<d) 
mentioned in the preceding paragraph. The evolution of the 
glassy dynamics of m-toluidine was measured for tempera- 
tures decreasing from 295K to 250K. The response functions 
could be fitted well by the scaling-law results predicted by the 
leading-order asymptotic formulas for the MCT bifurcation 
|0]. The analysis implies a critical temperature T c near 220K 
and a critical exponent a near 0.3. However, lowering the tem- 
perature to 225K, the critical power-law decay %(t ) °< l/t l+a 
was not observed |8]. The negative slope of the measured 
logx-versus-log? curve is not 1 + a, rather it is a number 
smaller than unity, say 1 — b' . Decay laws %{t) °c 1 Jt l ~ h with 
exponents b' around 0.2 have been identified by Cang et. al 
121 for a number of other van-der-Waals liquids. The OKE 
response of salol was studied by Hinze et. al, and the data 
were shown to be consistent with the known MCT scaling-law 
formulas for the temperature decreasing from 340K to 266K 
The glassy response for T = 257K was measured with 
an impressive accuracy for times increasing up to 0.5 /js; and 
the dynamics for t > lOOps displays the behavior expected 
from the leading-order asymptotic results of MCT. However, 
the glassy dynamics for 2ps < t < 20ps manifests itself by a 
%(t) oc \/t decay. The corresponding correlator shows a log- 
arithmic time dependence. All OKE response functions mea- 
sured so far demonstrate a glassy dynamics that cannot be de- 
scribed by the general leading-order asymptotic formulas for 
the MCT bifurcation in the regime t < 30ps. In a recent alter- 
native approach II ill is was possible to fit some of the OKE 
data for t > 30ps. 

It was argued recently that the new facets of glassy dynam- 
ics discovered by OKE spectroscopy l8,l9t Holl can be under- 
stood as generic implications of the Cole-Cole law for the crit- 
ical dynamics 111 211 . In the following, this statement shall be 
explained in detail. Section ITT1 summarizes the equations of 
motion and the known scaling-law results of MCT. It is ex- 
plained in Sec.|ni]that an asymptotic expansion of the modu- 
lus can have a much larger range of validity than the expan- 
sion for the susceptibility; the Cole-Cole law - introduced in 
1941 as empirical law lfl3tl - is derived in full generality from 
the microscopic equations of motion. The relevance of these 
results is demonstrated for a schematic model in Sec. HVI us- 
ing parameter values that describe the mentioned OKE data. 
The interplay between the Cole-Cole peak and the a-peak is 
investigated in Sec|V]where for specific parameter values the 
a-peak displays a wing. Section lyTl presents a conclusion. 



II. ESSENTIAL MCT FORMULAS 

A. Equations of Motion 

Within the basic version of MCT, the dynamics of the sys- 
tem is described by M correlators § q {t), q — 1, . . . ,M. These 
are real and even functions of the time t, which obey the initial 
conditions ty q (t = 0) = 1, dt§q{t = 0) = 0. The corresponding 
set of normalized response functions is given by 



(la) 



Laplace transforms map functions from the time domain, say 
F(t), in the frequency domain. They shall be used with 
the convention LJ[F(t)](z) = i Jq dt exp[izt]F (t) , z = tt> + 
/O. One gets LT[x q (t)]{z) = [1 + CO(j) c/ (co)], where ty q ((i>) = 
LT[(j) 9 (r)] (z). The normalized dynamical susceptibilities are 
given by 



X q ((0) = l+ax^cu). 



(lb) 



The loss spectra X q ((i>) = ^ m Xti{ & ) are related trivially to the 
fluctuation spectra (|>^(co) = Imcj)^©): Xq(<*>) = a><K,' (<*>)■ 

The Zwanzig-Mori formalism provides a fraction repre- 
sentation of §q{($) in terms of a fluctuating-force correla- 
tor M q (a): § q ((j)) = -l/{w-£^/[w+M ? ((u)]}. The pos- 
itive frequency H q quantifies the initial decay of the corre- 
lator § q {t) = 1 - [a,f] 2 / 2 + 0(? 3 ). Within MCT, a white- 
noise term v q is split off from the kernel M q ((o); the remain- 
der is represented in terms of a dimensionless function m q (t): 
M q {&) = iv q +Q^m q (u3); v q ^ 0. Here, m q (t) and m q (a) are 
related by Laplace-transformation. The fraction representa- 
tion is equivalent to the equations of motion 



d^ q {t)+v q d t § q (t) 



%(t)+ / dt'm q (t-t%4 q (t') 



0. 



(2a) 



The essential approximation in MCT is the expression of 
m q (t) as a polynomial J q of the correlators: 



m q{t) = ?q[$l(t),---AM{t)} 



(2b) 



There is no monomial contribution of order zero. The co- 
efficients of the polynomial are called mode-coupling coeffi- 
cients. They are the coupling constants of the theory and must 
not be negative. Within the microscopic theory, the polyno- 
mials are of second order and the coefficients are given by 
the equilibrium structure functions, which in turn are smooth 
functions of control parameters like the temperature T for the 
states considered. 

At the generic transition mentioned in the preceding sec- 
tion, the correlator's long-time limits depend singularly on 
the distance parameter e. For e < 0, fluctuations disappear 
for long times, (j) 9 (f — > °°) = 0. For e ^ 0, § q (t — > °°) = f 
< f q < 1. q — 1,...,M; the fluctuations arrest. Within 



3 



the microscopic version of MCT, the arrested part f q has the 
meaning of the Debye-Waller factor of the solid amorphous 
state. The leading-order variation with changes of £ for the 
arrested part is given by f q — f q °= ^/e, £ — * + , f q > 0, 
q = 1, . . . ,M. The regularity of the MCT equations implies the 
following. For £ tending to arbitrarily small values, there ap- 
pears an arbitrarily large time interval, where 9 (f) is arbitrar- 
ily close to f q . This critical arrested part f c q has the meaning 
of a plateau for the § q (t )-versus-logf curves for states near the 
transition. The corresponding plateau for the force correlators 
is m q (t -» oo) = /»' c = T c ijc , . . . j^] . At the transition, if q (co) 
and m c/ (co) exhibit poles —f q /(0 and — f q nc /{£>, respectively. 
Hence, there is a region of small |e| and small co, where G)+iv q 
can be neglected compared to £l q m q (co) . This region is the one 
for the MCT glassy dynamics. The fraction representation of 
the correlators simplifies to (^(co) = — l/[co — 1 /m £/ (co)]. This 
formula can also be noted as 



com 9 (co) = co<j) c/ (co)/ [1 + co<|> 9 (co)] . (3a) 

Equation Jlbi yields the equivalent expression for the dynam- 
ical susceptibility: 



X,(co) = l/[l-com,(co)] . (3b) 

Within the regime of glassy dynamics, [1 — cora^co)] has the 
meaning of a modulus for the response described by % q {t). 
The pair of Eqs. (I2bll3al can fix the solution only up to some 
overall time scale fo. The latter is determined by matching of 
the transient dynamics with the glassy dynamics. For further 
details and for a list of original papers, the reader can consult 
Ref. H. 



B. Scaling Laws 

There is a straight-forward recipe to calculate from the cou- 
pling coefficients at the transition point and from the critical 
arrested parts f q a number X, 1/2 ^ X < 1. It is called the ex- 
ponent parameter for the chosen transition point of the model 
under discussion. It fixes the critical exponent a, mentioned in 
Sec.|I] It fixes a further exponent b, < b ^ 1, which is called 
the von Schweidler exponent. The equation for the two expo- 
nents reads r(l -a) 2 /T{\ -2a) = X = T(l +bf/T{\ +2b). 
Parameter X also specifies a pair of equations for a pair of 
functions g±(i), which are defined for i > 0. The equations 
read ±1 + Xg±(i) 2 = (d jdt) / * di g±(i - i')g±(t'), and they 
have to be solved with the initial condition lirri£_^o t a g± (f ) = 1 . 
Up to corrections of order i", one gets for small t : 



g±(f«l) = l/f a . (4a) 

The function g+(i) approaches its long-time limit exponen- 
tially, g + (i ^> 1) = 1 /a/1 —X. The function g-(t) exhibits a 



power-law divergence for large i . Up to corrections of order 
1 /f h , one gets 

S_(f »l) = -,Bf fo . (4b) 

There are tables allowing the determination of a, b and B from 
a given X. There are also tables to determine g+(i) with an 
accuracy sufficient for all practical purposes 11511 . 

The functions g±(i) are the shape functions for the first 
scaling law of MCT. For E tending to zero, there appears a time 
interval of diverging length, within which A q (t ) = |(j) c/ (f ) —fjj\ 
is arbitrary small. The leading-order solution for the small 
parameter A q (t) yields the first scaling-law. The solution as- 
sumes the form 

Ut)=f«+h<iV\a\8±(t/tc), e^O. (5a) 

The amplitudes h q > 0, q = 1, . . . ,M are calculated from the 
mode-coupling coefficients at the critical point. The separa- 
tion parameter a is defined similarly: it is a smooth function 
of the control parameters and can be linearized close to the 
transition point, a = Ce + (e 2 ), with a positive coefficient C 
that depends on the chosen control parameter. The first critical 
time scale t a is given by the critical exponent a, the separation 
parameter a, and the time scale fo, which is defined by the 
short-time dynamics: 

f o = f /|a| 5 , 8 = l/(2a). (5b) 

The strong control-parameter dependence of the correlators 
near the plateau is described solely by that of the correlation 
scale and of the time scale f . 

From Eqs. J4al5al5bl . one gets the leading-order asymp- 
totic law for the decay of the critical correlator discussed in 
SecID § q (t) -f% = h q (t /t) a . From Eqs. (gbj |5al [5bJ one gets 
the von Schweidler law for the decay of the liquid correlators 
below the plateau: 

§q{t)=f q +h q {t/t & ) b , f a <f<? ', e<0. (6) 

Here, the second critical time scale of MCT reads f<y = 
f /[£ 1/fo |cTM, withy= l/(2o) + I /{2b). Up to errors of or- 
der |e|, the decay of the liquid correlators below the plateau is 
described by the second scaling law of MCT 

<M0=$?(*Ac), f a<f, £<0. (7) 

Here, § q (t ) is an £-independent shape function. Its initial part 
is given by von Schweidler's law: (f <C 1 ) = f q — h q t b . Cor- 
rections of order \/|ij modify the formula (0 for the below- 
plateau decay in a regime where § q {t) is close to the plateau. 
These corrections are described by Eq. l l5al for t ^ f c . The 
below-plateau decay is referred to as a-process, and Eq. Q 
formulates the superposition principle. Details of the deriva- 
tion of the cited results, references to the original work, and a 
comprehensive demonstration for the hard-sphere system can 
be found in Refs. Q d . 
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III. CRITICAL DYNAMICS AND COLE-COLE LAW 



B. Power-Law Solution for the Susceptibilities 



In this section, correlators and susceptibilities shall be dis- 
cussed for states at the transition point, say T = T C . Focusing 
on the range of validity, the essential formulas are analyzed 
for the correlators in lHI Al and for the susceptibilities in llll Bl 
The Cole-Cole law is derived in lHI CI 



A. Power-Law Solution for Correlation Functions 



The leading-order scaling-law formulas J4al5at yield for 
the critical correlators 



Ut)=fq+h q (t0/tY 



(8) 



A first question, whose answer is not implied by the results 
cited in Sec. Ill Bl concerns the range of validity of Eq. for 
short times. Let us define an onset time t * for the power law by 
the request, that Eq. l|8) describes § q (t) — fS within a relative 
error of, say, 10%, for t ^ t* For normal-liquid dynamics, 
one would expect the correlators to decay to zero for times 
around t m \ c . A second question to be discussed is: How can 
one describe the critical correlators in cases where t * exceeds 
fmic, i.e., when there is a gap between the end of the transient 
regime and the onset of the critical power law? 

The asymptotic solution l|8) can be extended to an asymp- 
totic series expansion in powers of t~ a . The result up to the 
next-to-leading term shall be noted as 



<M') = ./£+V'oA) fl [i+£?('o/O a ] 



(9a) 



Here, remainders which are given by (fo/f) 3 " times some 
power of ln(f/fo) are dropped. The Tauberian theorem yields 
an equivalent formula in the frequency domain: 



co(^(co) = -f q -T{\ -a)h q (-mt ) a 

-F(l-2a)h q K q {-mt ) 2a . 



(9b) 



In this formula, terms are dropped, which are proportional to 
times some power of In CO. The preceding two formu- 



co 



3a 



las are the starting point for the following derivations. It is 
a straight-forward procedure to calculate the correction am- 
plitude Kq from the mode-coupling coefficients at the critical 
point T = T c (3- From Eq. d9al . one can estimate an onset 
time of t */to — [10|a ? |]5. This is an estimate based on the as- 
sumption that higher-order expansion terms in Eq. d9al do not 
influence the results seriously for t ^ t*. Typically for many 
systems, the critical exponent a is around 0.3, and I /a > 3. 
Hence, t* depends sensitively on the correction amplitude K q . 
As a result, t * can vary considerably for different q. As a rele- 
vant example, let us cite the results for the density-fluctuation 
correlators of a system of hard spheres of diameter d 0: 
The time t\ referring to a wave-number g \d = 7.0 exceeds the 
time ?2 for qj_d = 10.6 by a factor of around 100. 



From Eqs. (fTall8t. one obtains the long-time result for the 
response functions up to leading-order corrections, 



X(f)fo = ah q (t /t) 



l+o 



(10) 



According to Eq. J9at . the leading order result in Eq. (|8ji is 
valid for t > t* but the onset time for Eq. dlOt is later, t ^ 
f*2'/", where typically 2 1 /" sa 10. Hence, the detection of 
the critical dynamics in its leading asymptotic form is more 
difficult for the response functions than for the correlators. 

From Eq. \9b\ , one arrives at the expansion formula for the 
absorptive part of the dynamical susceptibility 

4 (CO) = [r(l -a) S m(na/2)]h q ( m) y 

r » i (Ha) 

x [\+k a K q {m) a \ , 



k a = 2T(1 — a)cos(na/2) /X. 
The leading-order power-law result reads 



(lib) 



4(co) = [r(l - a) sin (jw/2)] h q {m ) a , (11c) 

which has an onset frequency of CO* = l/(t*kl/"); for CO ^ 
CO* the leading-order result in Eq. i 11 ci describes the critical 
spectrum with an error smaller than 10%. If A, decreases from 
1 to 1 /2, k a increases from 2 to near 5. For A, near a typical 
value of 0.7, k a is above 3, and the onset frequency for the 
CO a -law is about 30 times smaller than l/t*. Therefore, the 
detection of the critical power-law in the loss spectra is even 
more difficult then in the response function. 

Including the leading-correction terms in the asymptotic 
formulas, as noted in Eqs. l l9"al II lal i, is an obvious manner 
to extend the range of applicability of the analytic description 
of the dynamics. This is demonstrated comprehensively for 
the MCT for the hard-sphere system in Refs. 114 Il6ll . But 
there are cases, where this procedure does not lead to satisfac- 
tory results. The mean-squared displacement is an example, 
where the description of the increase towards the plateau by 
the analog of Eq. l l9"at cannot account for the glassy dynamics 



C. The Cole-Cole Law 

Equation ( I9bl i is an asymptotic expansion in terms of pow- 
ers of the small quantity t, = (— /cofo) fl that holds up to errors 
of Substitution of this expansion into Eq. Oat provides 
an analogous expansion for com f/ (co). Let us indicate the non- 
trivial parts of the coefficients by a superscript m; 



-C0m ? (co) = /; ,£ +r(l - a)h^(-mt ) a 

+ T{\-2a)h™K q n {-mto) 2a ■ 



(12) 
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Comparing coefficients of equal powers of i;, one finds: The correction amplitude reads K q c = K™(\ + f" ic )/ ), 

which is equivalent to 

f? c =f q /(l-®, f q =f? c /( l +f? C )- < 13a ) ^ q c =\+[{\-f q )/h q )k q /X, (16c) 

If the frequencies are so small that \K q c {— ;'co/co^)| <C 1, 
, n0 . t Eg. (1151 simplifies to the transparent expression of the Cole- 

\ =h q /{l-f q ) , h q = h q /(\+f q ) . (13b) Cole law, 



k'- = k cl + i[h q /{\~f q )], 
k q = k' q "-X[h'-/(i+f q ' c )]- 



(13c) 



Since the auto-correlation functions are normalized, f/ (f = 
0) = 1, one gets f q < 1. The kernel m q (t — 0) can have any 
positive value. Therefore, in principle, f' nc can be any posi- 
tive number. The different normalizations can be eliminated 
in the critical amplitudes by comparing the ratios h q J f q and 
h™/f' q nc : 



K/f? c ] = [h q /f q ]/(l-f q ). 



(14) 



These ratios determine the relative amplitude of the dynamics 
around the plateau within the first scaling-law regime. Equa- 
tion (EJ yields {if q {t)-f c q )/f c q = [V^VW*±('A<0. and 
a corresponding identity holds for m q {t). From Eq. dl4> one 
concludes that the relative amplitude of the first-scaling-law 
contribution is larger by 1/(1 — f q ) > 1 for the modulus than 
for the susceptibility; and 1/(1 — f q ) increases with f q , the a- 
peak contribution of the loss spectrum. Consequently, within 
the range of validity of the first scaling law, it is easier to mea- 
sure the first-scaling-law contribution for the modulus than 
that for the susceptibility. 

The expansion dl 2i can be substituted into Eq. j3bt in or- 
der to obtain an asymptotic expansion for the inverse of the 
critical susceptibility, i.e., for the modulus. The result shall be 
noted in the form 



x» = yfij i + Hco/O" + f% (-iq>/«£ 



, 2 a 



(15) 



The expression % q {(£>) is correct up to errors of the order 
C0 3fl . The three parameters specifying this formula can be 
expressed in terms of the coefficients in Eq. Mil and, via 
Eqs. H3i . in terms of the coefficients specifying the correla- 
tors in Eq. (|9aJ. One gets % c Q c q = 1/(1 +f' q nc ). This amplitude 
is the complement of the oc-peak strength f q of the normalized 
loss spectrum, 



- f 
J q 



(16a) 



The characteristic frequency entering the new formula for the 
susceptibility, is given by (co^o)" = (1 +f q c )/{h q T{\ - a)) 
or by 



a#b=[(i-yj)/(*,r(i- 



0)] 



L /a 



X q (v)=xZ q /[l + (-m/<*' q y] 



(17) 



The condition of validity for this formula means that the 
critical modulus can be described by the simple power 
law in the first line of Eq. dl 2i . The modulus spectrum 
obeys a formula analog to Eq. (II lc> : cora^(co) = T(l — 
fl)sin(7tfl/2)/i™(cofo) a . Using the 10% criterion from above, 
the last term in Eq. (I12> can be dropped for frequencies 

cof < [V (l0r(l -a)|£™|)] The correlator correspond- 
ing to the loss spectrum in Eq. dl 7i is given by the Mittag- 
Leffler function of index a: M u (x) =L,7=o x 'Vr(l +na), 



Ut)=f q +Xo q M a [-(m c q : 



(18) 



(16b) 



The Mittag-Leffler function can be calculated efficiently by 
Fourier-back transformation of %' q ((ti) / CO. 

Equation d!7t describes a very broad peak for the loss spec- 
trum x'^co) with a maximum located at co = co^. The spectrum 
is invariant under the interchange of (co/co£) to (co^/co): the 

(co)-versus-logco curve is symmetric. The width of the peak 
increases strongly with decreasing exponent a. For a 0.3, CO 
has to increase by more than a factor of 10 4 in order to scan 
the interval of frequency where % q ((d) %' q ((O q )/2. This sus- 
ceptibility formula known as Cole-Cole law, Eq. (I17> . was in- 
troduced as an empirical formula for dielectric loss spectra in 
glassy systems [13]. In particular, broad peaks located above 
the a peaks - often called P-peaks in this context - have been 
fitted by it. 

The Cole-Cole formula, Eq. (1171 . exhibits simple limits for 
small and large frequencies. For (co/co^)" <C 1, the Cole-Cole 
susceptibility reproduces the general power-law for the loss 
spectrum, Eq. dl lc> . % q ((o) « CO". Similarly, the correspond- 
ing Mittag-Leffler correlator reproduces the general long-time 
asymptote for the response, Eq. dlOI . For (co/co^) a 3> 1, the 
Cole-Cole spectrum describes a critical spectrum which de- 
creases with increasing frequency: % q (<£>) « l/co a . This is 
similar to a von Schweidler-law spectrum. The correspond- 
ing Mittag-Leffler correlator decreases according to the law 
(])(f) — const. °< —t". This yields a response function %(t) oc 
1 /t x with x = 1 — a < 1. Such behavior is consistent with the 
one detected by the OKE results for van-der Waals liquids 1 9] . 

The crossover of the critical spectrum from the low- 
frequency wing to the high-frequency wing of the Cole- 
Cole loss peak has a counterpart for the Mittag-Leffler cor- 
relator that is most easily seen in a semilogarithmic plot. 
For (t(O q ) a <§C 1, the [0(f) — const.] -versus-logf curve is bent 
downward, as known for the von Schweidler curves. For 
(fcojp" 3> 1, there is the critical power-law decay which shows 
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up as an upward-bent curve. Hence, the 0(f )-versus-logf 
curve exhibits an inflection point. For an exponent a around 
0.3, the curve is nearly straight for an increase of log 10 f by 
a factor of about 3. In this case, there is nearly-logarithmic 
decay of the critical correlator for a time variation over three 
orders of magnitude. 

It depends on the size of 1 / CO^ relative to the microscopic 
time scale f m ; c which part of the Mittag-Leffler correlator 
dominates the critical dynamics in the time range of inter- 
est. The crucial question concerns the range of validity of 
the leading-order result for the susceptibility and for the mod- 
ulus. These ranges are determined by the correction ampli- 
tudes. If \K q \ is large and \K'"\ small, the Cole-Cole for- 
mula has a larger range of applicability than the general sim- 
ple power-law formula (II lcl . But, if \K™\ is larger than \K q \, 
Eq. M lei is a better approximation for the critical decay than 
the Cole-Cole spectrum. Since the brackets in Eqs. Jl 31 are 
positive, there are two obvious results. If K q is positive, there 
holds K™ > K q . In this case, the leading-order result for the 
susceptibility is a better description of the critical spectrum 
than the Cole-Cole spectrum. If K' n is negative, there holds 
\K™\ < \K q \ and the application of the Cole-Cole law is su- 
perior to the one of the general leading-order result for the 
loss spectrum. For negative correction amplitudes K q , there 
is a trend for cancellation of the two contributions to K'" in 

H 

Eq. (II 3cl . It is a generic situation, that \K n q \ is smaller than 
\Kq\. Hence, for K q < 0, it is expected that the Cole-Cole law 
is a good description of the critical dynamics. 

A detailed discussion of the MCT equations has shown that 
a large arrested part f c q implies a large negative correction am- 
plitude K q for the critical decay J3- In this case, the onset 
time f * may exceed the microscopic time scale f m ; c by several 
orders of magnitude. The time range for the applicability of 
the leading-order asymptotic description of the critical decay 
may be outside the accessible range of the existing spectrome- 
ters. But, this is the situation, where \K' q \ may be so small that 
the Cole-Cole law or the Mittag-Leffler correlator can provide 
an adequate description of the critical dynamics. 



IV. APPLICATION TO SCHEMATIC MODELS 

The preceding results shall be demonstrated in a simple 
schematic model for two different cases where the model pa- 
rameters are adjusted to reproduce the data measured for ben- 
zophenone (BZP) and salol I^HoIIH 

The simplest MCT models deal with a single correlation 
function only. For these models, the formulas in the preced- 
ing sections simplify since the label q can be dropped. The 
correlator, the critical arrested part, and the critical amplitude 
shall be denoted by 0(f), f c , and h, respectively. The model 
is specified by the frequency Q. and the friction coefficient v 
in the equation of motion d2al . Furthermore, there are the 
non-negative coefficients Vi, I ^ 1, which specify the mode- 
coupling monomial of order / for J in Eq. ( I2bl . In applica- 
tions for data descriptions, the specified numbers are consid- 
ered as smooth functions of the physical control parameters. 



The simplest model for the polynomial J , which can repro- 
duce all possible values for the exponent parameter X, is given 
by the following formula for the fluctuating-force kernel 

m(f) = v 1 0(f) + v 2 0(f) 2 . (19) 

The two coupling constants vi and v 2 specify the state of the 
system by a point in the first quadrant of the vi-v 2 plane, cf. 
lower inset in Figs.^andll] The points (vj,v|) for generic 
fold bifurcations are located on a piece of a parabola. The 
position of the specific transition point can be characterized 
by the value for X 11711 : 

v\ = (21- l)/X 2 , v c 2 = 1/A 2 , f = l—X,h = X. (20a) 

The separation of some point (vi , v 2 ) from the transition point 
(vj,v|), see upper insets Figs.^ an dE] is given by 

a=[vi + v 2 (l-A,)]A,(l-A,), v l . 2 = v ia -v\ 2 . (20b) 

The arrested part of the correlator for the glass state reads 
0(r -> 00) = h[a/(l - X)} '/ 2 + O (a 3 / 2 ). The splitting of this 
value in an amplitude h and a remainder is not unique. The 
cited value for the critical amplitude h follows the conventions 
made in the preceding literature 11411 . The general expression 
for the correction amplitude K for M — 1 models [18] is spe- 
cialized easily to K = k(x), where the functions k(x) is defined 
by 

k(x) = -r(i -x) 3 /[xr(i -3x) -r(i -*)r(i -2*)] . (21) 

Function k(x) increases monotonically with increasing X: 
it decreases with increasing x. For X = 1/2, i.e., for a — 
0.395 . . ., one gets k(o) = —0.169 .... For X approaching 1, 
a tends to zero and k(a) diverges. For a = 1/3, X = 0.684 . . . , 
the correction amplitude vanishes. All higher-order correc- 
tions for the critical decay outside the transient regime van- 
ish as well for this special value of X fiUl . Therefore, for A, 
near 0.7, the simple power-law formulas of the leading-order 
asymptotic-expansion theory like Eqs. l|8] ^| II lcl describe 
the critical dynamics very well. 

For the specified model, there holds f c $C 1 /2. Hence, this 
model cannot be used to describe glassy-dynamics data for 
systems where the oc-peak loss spectra have a weight f° which 
exceeds 50% of the total weight. Moreover, the relation be- 
tween the exponent parameter X and the critical arrested part 
f c , as formulated by Eq. ( I20al . is an artifact of the model. The 
minimum requirement for a schematic model for data analy- 
sis is the freedom to adjust X and f c independently. This goal 
can be achieved by introducing a second correlator. The first 
correlator (j) 9= i (?) = 0(f) is used as a caricature of the density- 
fluctuation dynamics. It provides the exponent parameter. The 
second one describes the dynamics of some probing variable, 
say, A. The correlator shall be denoted by 9= 2 (f ) = 0a (t) ■ All 
quantities referring to this correlator shall be indicated by an 
index A. The equation of motion ( l2al is specified by Q.a > 0, 
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Va 0, and a kernel rriA(t). The two frequencies quantify the 
transient dynamics. The kernel is a polynomial ( 12bi of the 
two correlators involved. The simplest model describing the 
coupling of the probing variable to the density-fluctuation is 
given by 12011 



m A (t) =va<K0<I>a(0' 



(22) 



The coupling of the probing variable to the density fluctua- 
tions is quantified by va > 0. Since there is no influence of 
the dynamics of the probing variable on the density dynamics, 
the same scaling-law function g±(t/t a ), Eq. ( Bat , describes 
the leading-order near-plateau dynamics of 0a (0 as it does 
for that of 0(f). Notice in particular that the parameters Q.a, 
Va, and va do not modify the time scale to. In applications for 
data descriptions, the model parameters Q.a, Va, and va are 
to be considered as smooth functions of the physical control 
parameters, cf. lower insets in Figs.[0and|4] 

The loss spectrum % A (co) develops a (3-peak upon increas- 
ing the mode-coupling coefficient va [21]. For this spe- 
cial case and in the limit of infinite va, the Cole-Cole sus- 
ceptibility and the Mittag-Leffler correlator have been de- 
rived for the critical dynamics of the probing variable in Ref. 
llisil . The model specified by Eqs. ( 1191 I22> was used re- 
peatedly for the description of experimental data 
l23l l24l I25I |2(J. Large data sets for the evolution of the 
glassy dynamics of propylene carbonate have been analyzed 
in Ref. 1 27 ] ; results measured for different probing variables A 
obtained by neutron-scattering, depolarized-light-scattering, 
and dielectric-loss spectroscopy have been fitted by a com- 
mon first correlator 0(f). The different probes have been char- 
acterized by adjusting the parameters for the second correlator 
0a(O only. 

Up to order (fo/f) 2a , the kernel m.A(t) can be calculated 
from Eq. i22\ substituting the expression d9al for 0(f) and the 
analog expression for the second correlator, 



♦a(«) =fA + ha{to/t)" [1 +K A (t /t) a ] 



(23) 



Laplace transform yields cowm (co) as an asymptotic power se- 
ries in the small parameter i; = (— /cofo) a . The coefficients are 
linear functions of f jj, Iia , and Ka- This is substituted on the 
left-hand side of Eq. ( l3al for q = 2. Expanding the right-hand 
side up to orders one can compare the coefficients in order 
to arrive at: 



/a = 1-[1/(/ c va)], h A = X/{f 2 v A ). 



(24) 



The correction amplitude shall be noted in the convention 
of the general theory 0: 



K A =K(a)+K A , 



Ka — Xv a 



1 



va(I-X) 



-X 



(25a) 



(25b) 
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FIG. 1: OKE response functions measured for BZP for T/K = 
251, 260, 290, 320 Hi (full lines from bottom to top) and fits by the 
schematic-model functions %aM 1 12] (dotted lines). The straight 
dashed line has slope —0.80. The dash-dotted lines labeled sc are 
the approximations by the scaling functions Eq. <28t . the thin dashed 
lines are fits with Eq. 1281 for freely adjusted s c and t a . The insets 
show the separation parameter a (right) and the fitting parameters 
(left); two additional state points (+ and x) indicate extrapolations 
for which solutions are shown in Fig.0 



This leads to the parameters for the Cole-Cole susceptibility: 



1 c oa = 1/(/ c va) , co^f = [(1 -X)/(Xr(l -a))} 



1/" 



(26) 



Remarkably, the comparison of Eq. ( I26> with the general ex- 
pression in Eq. dl6b> shows, that the Cole-Cole frequency 
(0 A does not depend on the coupling coefficient va- CO^fu de- 
creases with increasing X from 0.37. . .for X = 1/2 to about 
10~ 4 for X near 0.88. For X around 0.7, co^f is about 0.02. 



A. Critical Relaxation for a Small Cole-Cole Frequency 

Figure[2reproduces OKE-response functions measured for 
benzophenone (BZP) |9] and fits to these data by the response 
%a(0 calculated for the model defined in the preceding sec- 
tion. The fit parameters for the mode-coupling coefficients 
v\, V2, and va are specified in the lower inset. The analysis 
shall be done by anticipating a transition point for X = 0.70, 
which implies the exponents a = 0.33 and b — 0.64. The cross 
in the inset is close to this bifurcation point. The o-versus- 
T diagram is shown in the upper inset; and extrapolation to 
= suggests the critical temperature T c = 235K with an es- 
timated uncertainty of ±5K. The 25 IK result exhibits the ex- 
pected von Schweidler decay, log%(f ) = const. — (1 — b)\ogt 
for times between about 0.3ns and about 6ns. Further details 
can be inferred from Ref. 1 12]. 




log 1Q t[ps] 

FIG. 2: Correlation functions for BZP. The lower panel exhibits the 
correlators ^(f) underlying the fits in Fig. Q (heavy full lines from 
left to right for T/K = 320, 290, 260, and 251, respectively). The 
light full line with label c shows the critical correlator calculated for 
X = 0.7 and v A = 30. The dotted line is the leading asymptotic law, 
Eq. JS}, with the time scale fo = 0.3775ps. The dashed line marked 
vS shows a von Schweidler law, Eq. J6j, with b = 0.64 and a time 
scale adjusted to match the T = 25 IK curve. The dashed line labeled 
A: is a fit by the Kohlrausch law, $ A (t) = /^exp[-(r/t)P], p = 0.91, 
with x adjusted to match the T = 260K curve. The times to and 
1/co^ = 15ps are indicated by arrows. The upper panel shows the 
critical decay (c) together with the leading-order (t~ a ) and next- to- 
leading-order (t~ 2a ) asymptotic power-law solution, Eq. <9at . The 
dotted curve labeled cc displays the leading-order Cole-Cole solu- 
tion, Eq. {TJQ, Wlt h Xo = 1/9- The points of 10% deviation of the 
approximations t ~ a , t ~ 2fl , and cc from the critical decay are indicated 
by the diamond, circle, and triangle, respectively. 



The MCT results describe the measured evolution of the 
glassy dynamics adequately for times exceeding f m ; c = lps 
with two exceptions. The T = 25 IK data exhibit a pronounced 
oscillation near 1 . lps, while the calculated curve shows this 
oscillation near 0.9ps. Furthermore, there are some deviations 
between the theoretical results and the data for the T = 260K 
results for times around 10ns. The four correlators 0a (0 and 
the four loss spectra yjlit), which are shown in the lower 
panels of Figs. [2] and [3] respectively, are the MCT results 
corresponding to the MCT responses JCa(0 shown in Fig. [2 
The data fits have been calculated for fixed relative values of 
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FIG. 3: Spectra for BZP. The lower panel shows the susceptibil- 
ity spectra for the correlators from the lower panel of Fig. [2] The 
leading-order critical spectrum, Eq. II lei , is the dotted straight line 
marked by co fl . A dashed straight line with slope —0.2 is fitted to the 
T = 25 IK spectrum to indicate an CO law for b' = —0.2. In the 
upper panel, the dashed curve marked co 2a shows the approximation 
by Eq. II lai . The approximation by the Cole-Cole function (cc), 
Eq. 1 171 . is shown dotted, the inclusion of the correction, Eq. I15i 
yields the dashed curve labeled ccc for co^ = 67ns~'. The points 
of 10% deviation of the approximations co~ a , co~ 2a , and cc from 
the critical spectrum are indicated by the diamond, circle, and trian- 
gle, respectively. The frequencies co^ and v = co/(2ji) = ITHz are 
marked by arrows. 



the four frequencies, which specify the transient dynamics: 
Q.A — £2, Va = V = 5£l. The scale Q. for the fits is chosen dif- 
ferent for different temperatures T 11211 . The curves in Fig. [2] 
and[3]are presented with a ^-independent scale corresponding 
to that one used for the fit for T = 25 IK: Q. = 1 .67ps~ 1 . With 
the mentioned reservations, these curves can be considered as 
the measured quantities for the OKE-probing variable A. The 
full lines marked c in Figs.|2]and|3]show the critical correlator 
(|>^(?) and the critical loss spectrum %'a ( <m )' respectively, for 
the transition point specified by X = 0.70 and va = 30. The 
correlator relaxes to the plateau f\ = 8/9 with a critical am- 
plitude h A = 7/27. 

The second scaling law for the decay below the plateau of 
the correlators, Eq. @, implies a corresponding scaling law 
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for the response: 



7U(t)=7U(t/t' a )/t' a 



(27) 



The control-parameter independent shape function reads 
%a(J) = — d§ A (i)/di, The reader can check this superposition 
principle for the results in Fig.[^as follows. The logXA(0~ 
versus-logf curve for T = 320K and t ^ 2ps can be translated 
so that it collapses with the curves for T /K = 290, 260, and 
251 for times exceeding 4ps, 90ps, and 290ps, respectively. 
Equivalently, the four loss spectra in Fig.[3]for % A (co) ^0.1 are 
connected by the superposition principle: % A ((ti) — % A ((£)t' a ) . 
The von Schweidler law describes a part of the high-frequency 
wing of the loss peak as shown by the straight dashed line 
with label vS for the T = 25 IK curve. The Kohlrausch law 
for a stretching exponent p = 0.91 describes the upper part of 
the oc-peak as shown by the dashed line with label K for the 
T = 260K result. 

There is an interval of times between the end of the tran- 
sient dynamics and the start of the von Schweidler decays 
that deals with the relaxation towards and through the plateau 
f A . The dynamics causes the loss spectra for % A ((o) ^0.1 and 
co/ (2k) < O.lTHz shown in Fig. [3] The leading-order asymp- 
totic solution of the MCT equations of motion deals with this 
part of the dynamics together with the von Schweidler-law 
part by Eq. d5al . This formula yields the first scaling law for 
the response 



%A{t)=h A s a x(t/t a ); 



(28) 



The shape function reads %(i) = —dg±(i)/di, T ^ T c . For 
T > T c , Eq. J4bi yields the von Schweidler response for large 
rescaled times t = t/t a : 1) = B ■ b/t x , x = 1 - b < 

1. For small rescaled times, Eq. J4al yields %(t < 1) = 
a/i x , x = 1 + a > 1. The crossover from one asymptote to 
the other occurs for times near f . The light dashed lines in 
Fig.[T]exhibit fits of the data by Eq. (|28j for X = 0.70. How- 
ever, both scales s a and t a are adjusted with the aim to achieve 
a good match in the von Schweidler-law regime. The l/t 1+a - 
law is not exhibited by the results in Fig. [2 For T = 320K, 
the crossover time f is located near lps, i.e., it is within or 
close to the transient regime. For the other three temperatures 
and times within the interval f m i c < t < t a , the measured re- 
sponse and the calculated functions % A (t) are below the values 
of h A Sai(t/t„). 

The straight dashed line in Fig.[2has a slope of —0.80. It 
demonstrates a pseudo-von Schweidler decay of the T = 25 IK 
results for 2ps < t < 20ps: § A (t) -const. °< -t b ' , b' =0.20. 
The 0a (f )-versus-log/ curves in Fig.|2]for t > 2ps approach 
and cross the plateau as downward-bent curves before enter- 
ing the von Schweidler decay regime. There are no inflec- 
tion points of the curves for some time near f a , as implied by 
Eq. (E3. 

The first scaling law describes a loss minimum of some 
value % m ; n at some position CO m i n . The shape of the logx^(co)- 
versus-logco curves is independent of the separation parame- 
ter o and the scales fix the position of the curve: Xmin x \/p[> 



COmin 1 /ta- The minimum is due to the crossover from the 
von Schweidler wing of the oc-peak, %"((£>) x l/co fo , to the 
power-law asymptote for the critical dynamics, x"(co) « co". 
For most practical purposes, the minimum can be approxi- 
mated by the interpolation formula [28]: 



XAM/Xmin= ^(C0/C0 m i n ) a +fl(C0 m i n /C0) / 'J / (a + b) . (29) 

However, the loss spectra in the lower panel of Fig, [3] do not 
exhibit loss minima, which can be described by the interpo- 
lation formula for the leading-order asymptotic result. The 
minima for frequencies CO near lps -1 do not depend sensi- 
tively on a. With decreasing \a\, the minimum position even 
shifts slightly upwards rather than downwards. With increas- 
ing frequency, the high-frequency part of the a-peak exhibits a 
crossover from the von Schweidler decay, %"(co) 1 / co'', to a 
pseudo von Schweidler decay, %"(03) x 1 /(O b . For T = 25 IK, 
this high-frequency wing extends from CO « 0.02ps~' to about 
0.4ps-' with b' = 0.2, as is demonstrated by the dashed 
straight line. The observed minimum is due to the crossover 
from the 1 / co /j wing to the spectral peak for the normal-liquid 
dynamics. The latter is located near v — ITHz. 

The scaling-law fits, which are shown in Fig.[Jby the light 
dashed lines, are misleading because the scales s a and t a have 
been adjusted freely. Calculating these scales from Eqs. J5bl 
gnB|E3|l, one obtains for T — 251K and 260K the dash-dotted 
lines marked by sc. The two distance parameters e = (T c — 
T)/T c are -0.07 and -0.11, respectively. These values are too 
large for a leading-order asymptotic formula to be applicable. 

The correlator 0a (f) for the lowest temperature under dis- 
cussion is close to the critical one for t up to about 30ps, as 
is shown in the lower panel of Fig. [2] Hence, the T — 25 IK 
curves exhibit critical glassy dynamics for the times between 
lps and 30ps. From Eqs. d25l one obtains the correction am- 
plitude Ka = —1.51. This large value yields an onset time 
t\ for the f~"-law which is beyond t — 10 3 ps as is marked 
by the diamond in the upper panel of the figure. Including 
the leading correction term yields the approximation by the 
dashed line denoted ?~ 2a . It improves the description of the 
critical decay so that it can be understood for times larger than 
about t = 30ps. But the expansion formula fl23i cannot be used 
to describe that part of the critical decay, which is measured 
for BZP and described by the solution of the two-component 
schematic MCT model. 

Substituting the cited values for A,, f c A , and h A into Eqs. dl3al 
I13bt . one gets the plateau for the modulus, f™ c = 8, and the 
critical amplitude h™ = 21. Notice that the relative strength of 
the first scaling law amplitude for the modulus is much larger 
than that for the correlator: h' A /f A " c — 2.6 versus h A /f A = 
0.29. From Eq. dl3cK one obtains the correction amplitude 
K A = 0.124. As expected from the discussions of Sec.|^ the 
onset time t' A * for the leading-order formula for the modulus, 



m A (t)=fr + h A (t /t) a 



(30a) 



is smaller than t A by more than a factor 1000. As a result, see 
the upper panel of Fig. [2] the equivalent formula Jl 8I > for the 
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correlator, 



M*)=/A + (l-/!)M a [- (m A ) a 



(30b) 



describes the probing variable A for all times exceeding t — 
lps. The leading-order asymptotic result for the modulus ex- 
plains the response for T = 25 IK quantitatively within the 
interval lps^ t ^ 30ps. 

Figure shows that the leading-asymptotic description of 
the critical loss spectrum, Eq. (II lei , would be relevant for the 
explanation of the loss minimum only in cases with co m i n < 
2 • 10~ 5 ps~'. The dashed line labeled co 2a exhibits the asymp- 
totic expansion for the critical loss up to the leading correc- 
tion, Eq. i ll 1 at . This formula is relevant for co < 10~ 3 ps _1 . 
Even this expression is unsatisfactory for the discussion of the 
BZP results because of the large correction amplitude. How- 
ever, the Cole-Cole spectrum describes the critical loss rea- 
sonably for CO < 0.2ps _1 . The correction amplitude for the 
Cole-Cole law, Eq. iflgcl is very small, K C A C = 0.076. The 
Cole-Cole law with leading correction, Eq. dl5> . is shown 
in the upper panel of Fig. |5]by the curve marked ccc, and it 
slightly above the leading-order result in the regime of large 
frequencies. 

For the system under study, the Cole-Cole frequency reads 
C0 A = 0.067ps _1 . This frequency CO A /(27T,) » 10GHz is small 
compared to the ITHz scale for the normal-liquid dynam- 
ics. The value — logco A is in the center of the logf interval 
studied by the results of Fig. ^ Therefore, the critical spec- 
trum relevant for the understanding of the data differs quali- 
tatively from the leading-order power-law formula. As a re- 
sult, the loss minima in Fig. [3] which are caused by the su- 
perposition of the critical spectra for CO near and above co A 
and the von Schweidler-law wing of the a-peak, differ dras- 
tically from the general low-frequency shape described by 
Eq. \29\ . For T = 25 IK, the superposition yields the pseudo- 
von Schweidler law manifested as a wing specified by the 
exponent b' = 0.20. This wing will be analyzed further in 
Sec. El 



B. Critical Relaxation for an Intermediate Cole-Cole 
Frequency 

Five OKE response functions measured for salol | loj] and 
fits by the functions % A (t) of the above explained schematic 
model are reproduced in Fig. [4] The filled dots in the insets 
specify the mode-coupling coefficients v\, v%, and v A used for 
the calculations. Further details can be found in Ref. 0, 
where, however, the value for Q. was misprinted. The se- 
ries of states (vi,V2) extrapolates to a transition point with 
X = 0.73, which implies a critical exponent a — 0.31 and a 
von Schweidler exponent b = 0.59. The extrapolation of the 
O-versus-T parameters to a = suggests a critical temperature 
T c = 245K with an estimated uncertainty of ±3K. The MCT 
correlators 0a (f) an d l° ss spectra % A (03), which are equivalent 
to the response Xa(0 in Fig.[4l are shown in the lower pan- 
els of Figs.|5]and|6] respectively. They can be considered as 
the measured results for the glassy dynamics of salol, since 
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FIG. 4: OKE-response functions measured for salol for T/K = 
247, 257, 270, 300, 340 (full lines from bottom to top) (ujl. The 
dotted lines are fits by the schematic-model response %a(0 calcu- 
lated with a = 2Q A = 10v A = 15.9PS- 1 [Tj. Symbols and curve 
styles are the same as in Fig.Q The slope of the straight dashed line 
is -1.15. 



the calculated and measured curves in Fig. |4] agree outside 
the transient regime. This holds with the reservation, that fits 
and measurements for the T = 270K results exhibit some dis- 
crepancies for times exceeding 20ns. The lines with label c in 
Figs.|5]and[6]show critical correlators and critical loss spectra, 
respectively, calculated for A, = 0.73 and v A = 55. The corre- 
lators relax to the plateau f A — 0.93 with a critical amplitude 
h A = 0.18. 

The test of the superposition laws for the long-time relax- 
ation parts of the correlators is left to the reader. This sec- 
ond scaling law of MCT relates the a-relaxation peaks for 
the loss spectra for % A (03) ^0.1. Kohlrausch-law fits and 
von Schweidler asymptotes have been added to the data for 
T = 257K and T = 247K, respectively, in order to empha- 
size that the results for the evolution of the below-plateau- 
decay process follows the familiar pattern. Notice that the 
response for T = 257K exhibits the von Schweidler-law de- 
cay, log% A (t) = const. — (1 — b)logt for the large time inter- 
val 0.1ns ^ t ^ 10ns. The measurement for this temperature 
permits a rather precise determination of the exponent b and, 
thereby, of X. 

The light-dashed lines in Fig. @] are fits of the data by the 
first-scaling-law expression J28t i with freely adjusted scales s a 
and t a . For T ^ 270K, these fits describe the data outside the 
transient regime up to times within the von Schweidler-law 
region. A more detailed documentation of this result can be 
found in Ref. 11011 . where also consistency of the fitted scales 
s a , t a , and t' G with the MCT-power-law formulas was demon- 
strated. The scaling-law results calculated with the MCT val- 
ues for the scales h A , s a , and t a in Eq. J28t i are shown in 
Fig- |H as dash-dotted lines marked sc for the T = 257K and 
T = 247K data. For T = 257K, e = -0.05 is so large, that 
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FIG. 5: Correlation functions for salol. The lower panel shows 
the correlators used for the response fits inFig.Elfor T/K = 
247, 257, 270, 300, and 340 (heavy full lines from right to left). The 
light full line shows the critical correlator calculated for X = 0.73, 
and Va = 55. The dashed line marked vS exhibits a von Schweidler 
law, Eq. for a time scale chosen to match the T = 247K cor- 
relator. The dashed line marked K is a Kohlrausch-law fit for the 
T = 257K curve with exponent (3 = 0.95. The dotted line exhibits 
the formula J8} with to = 0.00246ps. The upper panel shows the crit- 
ical decay (c) together with the leading- (t~ a ) and next-to-leading- 
order asymptotic solution <9a> (t~ 2a ). The dotted curve labeled re- 
displays the leading-order solution from Eq. i 1 8i with %q c = 0.067 
and co^ = 7.25ps. The points of 10% deviation of the approxima- 
tions t~ a , f~ 2 ", and cc from the critical decay are indicated by the 
diamond, circle, and triangle, respectively. The times to and 1 /(if A 
are marked by arrows. 



FIG. 6: Spectra for salol. The lower panel shows the susceptibility 
spectra %^ (co) for the correlators shown in the lower panel of Fig. [5] 
The leading-order approximation to the critical spectrum, Eq. il lal . 
is shown as dotted straight line marked co". The dashed line labeled 
sc' exhibits the interpolation formula <29t with the von Schweidler 
exponent b = 0.59 and the critical exponent a replaced by a" = 0.24. 
The upper panel exhibits the critical spectrum as full line with label 
c, the leading-order asymptotic approximation, Eq. J 1 lcl . (dotted, 
labeled co") and Eq. il lal . (dashed, labeled co 2fl ). The approximation 
by the Cole-Cole function for co^ = 7.25ps (cc), Eq. <17> . is shown 
dotted, the inclusion of the correction, Eq. <15t . yields the dashed 
curve labeled ccc. The points of 10% deviation of the approximations 
co - ", co -2 ", and cc from the critical spectrum are indicated by the 
diamond, circle, and triangle, respectively. The frequencies co^ and 
Co/(2ti) = V = ITHz are marked by arrows. 



Eq. (12 8 1 cannot account quantitatively for the scales ruling 
the von Schweidler-law decay. For T 270K the distance 
parameter e exceeds 10%; the leading asymptotic expansion 
formula for the plateau-crossing process is not applicable for 
such large distances of the control parameter from the critical 
value. In contrast, the scaling law ( 12 8 1 describes the T = 247K 
result well for the large time interval 0.04ns < t < 10ns. In 
this case, the distance parameter e = —0.008 is so small that 
the leading-order expansion result for the plateau crossing, 
Eq. ( Bal . accounts for a major part of the long-time response 
- readjusting the scales (/i^Sc), and f would not improve the 
fit. 



Figure [5] shows that the T = 247K response is close to 
the critical one for t ^ Ins. The straight dashed line with 
slope —1.15, which is shown in Fig.|4] demonstrates that the 
critical response exhibits a power-law decay Xa(0 x 1/f 1+ " , 
a' = 0.15, for 2ps < t < lOOps. This is equivalent to a power- 
law decay of the correlator: <|>a(0 —f% 05 l/ r " ■ I n qualita- 
tive agreement with the prediction by Eq. ( Bal . the crossover 
from the f~" decay to the — t b decay causes an inflection 
point for the (j)^(f)-versus-log/ curve in Fig. [5] The cor- 
responding crossover from the GT b wing of the a-peak to 
some co" spectrum causes the loss minimum located near 
log co = —2.5 in Fig. [6] The dashed curve marked sc 1 exhibits 
a description of this minimum by the function (co) / Xmin = 
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b((o/(Q m i n ) a +a"(cOmin/co) fe /(a" + b) which is suggested 

by Eq. (129b . For salol, as opposed to BZP, the minimum 
position decreases strongly with decreasing |e|, in qualitative 
agreement with the first-scaling-law results. 

The inconsistent exponents for T = 247K in the regime 
2ps < t < lOOps indicate that the leading-order asymptotic de- 
scription is still not applicable; the exponents a' = 0.15 ob- 
tained from the fit of the OKE data and a" = 0.24 obtained 
from the fit of the corresponding minimum are different and 
both are smaller than the correct value of a = 0.3 1 . The devia- 
tion of the critical correlator from the leading-order power-law 
result (|8j and dlOt within the specified time interval is caused 
by the large value for the correction amplitude: = —1.78. 
This value implies an onset time for the t~ a decay of be- 
yond 20ps as shown by the diamond in Fig. [5] The high- 
frequency part of the loss minimum is located in the region 
10ns _1 < CO < lps -1 . Figure [6] demonstrates that the co" law 
does not describe the critical loss there. Inclusion of the lead- 
ing correction terms for the analytic description of the critical 
dynamics, i.e., using Eq. ( l9at for the correlator and Eq. i ll 1 at 
for the loss spectrum, explains the result for t 0.6ps and 
CO ^ 25ns~', as demonstrated in Figs.|5]and|5J respectively. 

From Eq. ( I13ct . one derives the correction amplitude for 



V. COLE-COLE WING 



the modulus K' A " = 0.177. The leading-order formula OOal 
for the relaxation kernel describes the kernel niA outside the 
transient regime. The corresponding expression ( I30bt for the 
correlator and the equivalent Cole-Cole formula for the loss 
spectrum account for the critical glassy dynamics, as shown in 
Figs. |3 and [6] The Cole-Cole formula with correction term of 
amplitude K% =0.09 yields a slight improvement compared 
to the equation based on K L A = 0, as shown in Fig.|6]by the 
curve with label ccc. Equation (I16bt leads to the Cole-Cole 
frequency for salol: co A = 7.25ps~'. This value is close to the 
loss peak for the normal-liquid dynamics: co A f m ; c ~ 1. The 
part of the glassy critical loss spectrum which is relevant for 
the explanation of the data in Fig. |4] deals with the regime 
10~ 4 < co/co^ < 10 _1 . Within this frequency interval, the 
Cole-Cole spectrum increases smoothly with CO. Therefore, 
the leading-order formulas d 1 Ob or il Id describe the dynam- 
ics qualitatively. However, the true critical spectrum differs 
from its low-frequency asymptote, Eq. il lcl . Hence, the de- 
scription of the dynamics by the scaling law is not correct 
quantitatively. Within the specified interval, the logx^(co)- 
versus-log co spectrum can be approximated reasonably by a 
straight line, log (to) ps const. + a" log CO, a" = 0.24. As a 
result, the critical dynamics is approximated well by a power 
law decay specified by an exponent a" smaller than the critical 
exponent a. 

In the preceding section, the critical spectra have been dis- 
cussed for both small and intermediate Cole-Cole frequencies 
CO^; the case of a large Cole-Cole frequency is found in the 
depolarized-light-scattering spectra for CKN I3l I2911 . In the 
fits with a schematic model by Alba-Simionesco and cowork- 
ers I3(ll3 ill one can identify a Cole-Cole peak for the critical 
spectra, but the description with the power-law solution is su- 
perior if terms up to order CO 2 " are considered. 



Figures^|3] anc l|4H6ldemonstrate the scenarios for the evo- 
lution of the glassy dynamics for of <C t^ c and cojj w f mi ' c , 
respectively, for the probing variable described by the second 
correlator of a two-component schematic MCT model. Solu- 
tions for this model, which exemplify a scenario as shown 
in Fig. |3j have been discussed before by Cummins l32ll . 
He pointed out that the evolution of the wing phenomenon 
obtained from the model is similar to the one known for 
dielectric-loss spectroscopy for glycerol and several van-der- 
Waals liquids |33j]. He emphasizes also that the cited mea- 
surements refer to temperatures T below T c , while the cal- 
culations are done for T > T c . Furthermore, the experimen- 
tal data demonstrate the wing only for frequencies below 
10GHz« lO- 2 ^. So far, no wing has been reported which 
occurs in the two-decade frequency window adjacent to the 
microscopic excitation region. 

In order to describe the wing phenomenon quantitatively, 
one can extend the Cole-Cole formula so that the spectrum 
can be described also for small but non-vanishing separa- 
tion parameters a — Ce°^ (T c — T)/T c . The leading-order ex- 
pression for the plateau-crossing process of the fluctuating- 
force correlators reads in analogy to Eq. ( Bat : m q (t) = f" lc + 
h'qg±(t /t c ), e ^ 0. Laplace transformation yields COm 3 (to) = 
—f™ c + h'qC((x)). Function C(co) obeys a scaling law, 



C(co) = vNc± (cofo) , e ^ , 



(31) 



with the control-parameter independent shape functions 
c± (to) given by the Laplace transforms g± (to) of the shape 
functions g±(t): c±(co) = cog±(co). Substitution of com ? (co) 
into Eq. d3bi yields the desired result. Expansion of the right- 
hand side in terms of the small parameter C(to) and compar- 
ison with the result for % c/ ((0) following from Eq. (Bat repro- 
duces the relations ( I13all3bf for the critical arrested parts and 
critical amplitudes. One gets 

Z 9 (to) = (1 -yj) /{l - [V(l -f q )] C(to)} . (32a) 

Equation ( 13 2 at was found for the schematic model in Ref. lH8ll 
for the special case of the second correlator 0a (0 an d m the 
limit of infinite coupling strength va- It is shown here, that 
Eq. ( 13 2 at can be derived without such restrictions. The lead- 
ing correction to Eq. ( I32at can be calculated from the asymp- 
totic expansions in 1 14l ll6ll in a straightforward manner; how- 
ever, the evaluation of the result for fitting data is rather in- 
volved. Focusing on the critical decay, one can extend the 
preceding formula by adding the large-frequency part of the 
correction: 



x» = (i-4)/{i-[V(W£)]c(to) 

+^ e (-/coM / ) 2a }. 



(32b) 



This expression describes the small-e dynamics for C0f m i c < 
1 and frequencies extending down to the beginning of the 
von Schweidler-law decay. 
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FIG. 7: Asymptotic description of the wing in BZP. The heavy full 
line reproduces the BZP loss spectrum for T = 25 IK from Fig. [3] 
and the dashed one exhibits the spectrum of the extended Cole-Cole 
susceptibility, Eq. <32at . The triangles mark the points of 10% de- 
viation between these spectra. The dotted line sc is the first-scaling- 
law result for the susceptibility, and the line cc shows the Cole-Cole 
spectrum for the critical point. The dash-dotted line represents the 
second-scaling-law result for the loss peak. The thin full lines are 
solutions corresponding to the state points indicated by + and x in 
Fig.Q 



BZP exemplifies the case of such a small K c q c that the 
leading-order result for the modulus, Eq. d32ai . can be used. 
Figure [7] demonstrates that this result accounts for the T — 
25 IK spectrum for the large dynamical range 10~ 3 ps~ 1 < C0< 
lps -1 . The scaling-law result for the loss spectrum, X^(co) = 
/!aC^'(co), is shown as dotted line sc; it can describe only a 
part of the von Schweidler-law spectrum for co < 10~ 3 ps _1 , 
but is inadequate for larger frequencies. This observation for 
the dynamics in the frequency domain is equivalent to the 
one demonstrated in Fig. [2 for the results in the time do- 
main. The second-scaling-law result for the loss spectrum 
is exhibited as dash-dotted line. It follows from Eq. 0, 
X^(co) = (cof^)^(cof^), and accounts for the loss peak for 
co < 10~ 3 ps _1 . Combining the result for the second scaling 
law for the modulus with the result for the first scaling law for 
the loss spectrum explains the BZP spectrum for T = 25 IK 
for co < lps~' including in particular the a-peak wing. The 
small but systematic discrepancies between the loss spectrum 
in the structural-relaxation regime and the asymptotic MCT 
expressions result from the fact that the distance parameter 
£ = (T — T c ) /T c k, 0.07 is so large, that the asymptotic results 
for the scales still have noticeable errors while the shape func- 
tions already agree well with the nontrivial spectra. 

Starting from the state that fits the 25 IK data for BZP, a 
number of extrapolations are possible within the schematic 
model. Keeping all parameters but (vi,v'2) fixed, two addi- 
tional states shall be considered; they are indicated by x and 
+ in the left inset of Fig. ^ an d their spectra are marked 
accordingly in Fig. Q The first state point (x) is close 
to the transition point and is characterized by a = —0.003 



(e = —0.01). Solutions at this point are shown in Ref. lll2ll 
for %(t) and Xa(0- The respective spectrum %^(co) in Fig.0 
exhibits nearly-constant loss for 0.4 10~ 3 ps _1 < CO < 0.4ps _1 
within a 10% margin. In this frequency window the maximum 
around co^ indicates an emerging Cole-Cole peak; the mini- 
mum at co = 10~ 3 ps -1 is caused by the crossover of the high- 
frequency wing of the a-peak and the low-frequency wing of 
the Cole-Cole peak. This minimum is ruled by the first scaling 
law and the divergent time scale t a . The second state (+) is an 
interpolation between x and the solution for BZP; the separa- 
tion parameter is a = —0.01 (£ = —0.03), and the point is cho- 
sen similar to the ones used in Fig. 4 of Ref. |32]. The spec- 
trum in Fig.Qshows a wing, XaI 03 ) ^ (a "' 1 > f° r almost three 
orders of magnitude in frequency. Hence, the crossover be- 
tween a- and Cole-Cole peaks can be interpreted as a power- 
law wing CO -6 for some range in frequencies and control pa- 
rameters without fine-tuning. The evolution of the crossover 
progresses from a wing with gradually lower exponents b' to 
nearly-constant loss and eventually the emergence of a Cole- 
Cole peak. While the extrapolations presented in this para- 
graph are not at all implied by the fit of the data, the scenarios 
seem nevertheless possible for realistic parameter values. 

I I I I I I I 
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FIG. 8: Asymptotic description of the OKE data for BZP. The full 
line reproduces the BZP data for T = 251K from |9]; the dotted line 
reproduces the fit from Fig. Q The dashed curve, the dash-dotted 
curve sc, and the dotted curve cc respectively show the Fourier back- 
transforms of the extended Cole-Cole susceptibility in Eq. <32at . of 
the first scaling law, and of the Cole-Cole law in Eq. 1171 . cf. Fig.0 

To conclude the discussion of BZP, let us return to the OKE 
data for T = 25 IK as shown in Fig.|S] In addition to data and 
fit for this state, the Fourier backtransforms of the asymptotic 
curves from Fig. are displayed. While the asymptotic so- 
lution at the critical point (labeled cc) covers only less than 
a decade, the extended Cole-Cole susceptibility in Eq. J32ai 
describes the data over three orders of magnitude in time and 
thereby explains the t~° ,8 -law found empirically. The remain- 
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ing discrepancies are again due to the relatively large distance 
parameter for this state; Fig. 3 of Ref. 11211 shows prefect 
agreement with the Cole-Cole law for states closer to the tran- 
sition including state x from Fig. [2 
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FIG. 9: Asymptotic description of the minimum in salol. The full 
line reproduces the salol loss spectrum for T = 247K from Fig. |6| 
and the dashed one exhibits the spectrum of the extended Cole-Cole 
susceptibility, Eq. <32at . The triangles mark the points of 10% de- 
viation between these spectra. The dotted line cc is the Cole-Cole 
spectrum for the critical point, and the line sc shows the first-scaling- 
law result for the susceptibility; the diamonds mark the points where 
sc deviates by 10% from the loss spectrum. The dash-dotted line 
represents the second-scaling-law result for the loss peak. 

Different from the case for BZP, the lowest temperature 
available for salol, T = 247K, is described by a much smaller 
distance parameter, e « 0.01, and therefore shares a larger fre- 
quency regime, say CO < O.lps -1 , with the critical loss spec- 
trum, see the lower panel of Fig. |6] Being so close to the 
critical point, the second scaling law in Fig. [9]- labeled a - 
is indistinguishable from the loss spectrum for CO ^ 10~ 4 ps -1 . 
The first scaling law sc describes the solution up to around 
co < 0.005ps~' which is the same as for the solution at the 
critical point, cf. upper panel of Fig. |6] As for the critical 
spectrum, the leading-order approximation in Eq. (|33aJ im- 
proves the description to larger frequencies by two decades 
up to CO ss O.lps -1 before crossing over to the Cole-Cole peak 
which in this case is hidden under the microscopic excitations. 



VI. CONCLUSION 

Within the regime of glassy dynamics, susceptibilities 
%q((o) and moduli COm ? (co) are related by Eq. Obi : and at the 
critical point, susceptibilities and moduli are both described 
by the universal power law co". The range of validity of this 
power law is given by amplitudes for the correction CO 2 " - K q 
for the susceptibilities and Kl" for the moduli, which are re- 
lated by Eq. ( I13ci . For some parameter regions a large value 
of \K g \ can render the power-law expansion for the suscep- 
tibilities irrelevant, while the correction amplitude \K'"\ is so 



small that the universal power-law can be used for the descrip- 
tion of the modulus successfully. Typically, this occurs if K q 
is large and negative, and this can be expected if the plateau 
value ft; is high, cf. Eq. J13ci . In this case, the susceptibility 
can be described well by the expansion of its inverse, Eq. J15I . 
For vanishing correction K™ = K^ c — 0, formula, the suscep- 
tibility at the critical point is given by the Cole-Cole law 1171 . 
a formula first introduced in 1941 II 311 . 

The Cole-Cole frequency co^ in Eq. ( 1171 introduces a char- 
acteristic scale for assorting the critical dynamics into three 
categories. (1) If COS is large compared to the scale t m ^ c for 
the band of normal liquid excitations, the Cole-Cole suscepti- 
bility reduces to the leading-order scaling-law formula which 
is implied by Eq. (BaV There is a control-parameter sensitive 
minimum as discussed for Eq. ( 1291 that interpolates between 
the von Schweidler-law tail of the a-peak, « 1 / co 6 , and 

the critical spectrum %q((o) x co", where a and b are related via 
the exponent parameter X. As an example for such a scenario 
one can cite the molten salt CKN. While the leading-order 
result alone is not sufficient to describe the measured data 
quantitatively, the power-law expansion in Eq. (II 1 ai yields a 
satisfactory description that is superior to the Cole-Cole solu- 
tion ll29ll . (2) If (if is close to f^V, one encounters a scenario 
demonstrated in Fig.|9]for salol. There is a control-parameter 
sensitive loss minimum. It originates from the crossover be- 
tween the von Schweidler-law tail and the critical spectrum; 
and the critical spectrum approaches the maximum of the un- 
derlying Cole-Cole peak. A description by Eq. ( 1291 is possible 
only, if the exponent a is replaced by some effective one a" , 
which is smaller than a, cf. Fig. [6] (3) If C0^ is smaller than 
the microscopic time scale, say C0^? m j c ~ 0.05, we obtain the 
scenario seen in Fig.0for benzophenone (BZP). In this case, 
the von Schweidler-law part of the a-peak, %'q x 1 /to' 5 , crosses 

over to some flatter wing, « 1 /co fo , a < b' < b. This wing 
is caused by approaching the Cole-Cole spectrum for frequen- 
cies around the maximum at (i£. For yet higher frequencies, 
one encounters a crossover - from the high-frequency part of 
the Cole-Cole spectrum for structural relaxation to the spec- 
trum due to normal-liquid excitations - producing a different 
kind of a minimum. The position of this minimum is control- 
parameter insensitive. For the understanding of this minimum, 
both the scaling-law i2$l and the interpolation formula ( I29t 
are irrelevant. On the other hand, as shown in Sec. IV! and 
especially in Fig.0 it is possible to describe the wing rather 
accurately by the leading-order formula ( I32al >. If the a-peak 
is shifted to yet lower frequencies, the wing can give way to a 
separate Cole-Cole peak. 

All figures in this work have been prepared with model 
parameters that reproduce the OKE-response functions of 
benzophenone (BZP) and salol. Hence, Fig. implies that 
the BZP spectra for temperatures near 250K, when mea- 
sured by depolarized light-scattering in backward direction, 
should exhibit an a-peak wing for frequencies v between 
about 1GHz and about 100GHz. The crossover from the 
von Schweidler-law wing to the wing induced by the Cole- 
Cole peak is expected to occur around 1GHz. This crossover 
position depends sensitively on the temperature, because the 
von Schweidler-law relaxation depends sensitively on T. 
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The Cole-Cole law with correction, Eq. (I15> . has been de- 
rived in its fully general microscopic version in Sec.|inj only 
in Sees. IIVI and IV1 the theory was specialized to schematic 
models. While for the schematic models a number of param- 
eters can be fixed to describe experimental data, for micro- 
scopic models the static structure of the model system deter- 
mines these parameters uniquely. A first example for Cole- 
Cole dynamics in a microscopic model has been discussed re- 
cently for the mean-squared displacement Sr 2 (f) of the hard- 
sphere system 13411 : including the Mittag-Leffler function for 
the critical relaxation allows for an analytic description of the 
mean-squared displacement for the full range of the dynamics 



similar to BZP in Fig. In addition, the Cole-Cole dynam- 
ics was identified in the data measured by van Megen et. al 
l35ll where Eq. d 1 31 accounts for the data for an interval in 
time of three orders of magnitude adjacent to the transient dy- 
namics. Hence, the results discussed above seem relevant for 
both molecular and colloidal glasses and can be expected to 
facilitate more detailed investigations. 
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